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Bave  spectra  calculated  using  the  Paraaeterized 
nonlinear  Wave  Solution  developed  by  Le  Mehaute  et  al. 
(1984)  are  coapared  with  field  data  acquired  at  Leadbetter 
beach t  Santa  Barbara ,  California.  The  paraaeterized  solu¬ 
tion  satisfies  the  nonlinear  free  surface  boundary  condi¬ 
tions  to  a  specified  degree  of  accuracy  and  is  expressed  in 
terns  of  a  converging  truncated  Fourier  series.  The  wave- 
number,  surface  profile  and  nave  orbital  velocities  are 
deterained  by  the  vave  height  and  vave  period  at  the  local 
depth  cf  water.  Spectral  coaponents  are  coapared  between 
the  aodel  results  and  field  data.  Good  agreeaent  is  observed 
for  waves  corresponding  to  Ursell  nuabers  ranging  froa  25  to 
75.  For  large  Ursell  nuabers  (strong  nonlinear  effects)  the 
paraaeterized  model  underestiaates  the  data. 
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The  random  nature  of  ocean  waves  is  most  .often  described 
using  the  linear  spectral  model,  and  many  engineering  prob¬ 
lems  related  to  offshore  and  coastal  structures  are  solved 
by  employing  the  concept  of  directional  wave  spectra,  fiaves 
in  the  sea,  however,  can  exhibit  nonlinear  properties.  A 
departure  of  surface  elevation  from  the  Gaussian  distribu¬ 
tion  with  positive  skewness  is  an  example.  Phase  velocities 
of  high  frequency  component  waves  not  satisfying  the  linear 
dispersion  relation  is  another  example.  Breaking  of  waves 
is  a  spectacular  example  of  nonlinear  behaviour  of  water 
waves. 

The  study  of  nonlinear,  water  waves  traveling  over  a 
horizontal  bottom  dates  back  to  Stokes (1847) ,  who  solved 
this  problem  in  the  form  of  a  power  series  in  terms  of  a 
small  parameter  '  6  *  related  to  the  average  wave  slope  ( 6  = 
ak,  where  a  is  the  wave  amplitude  and  k  is  the  wavenumber) . 
It  was  found,  however,  that  the  Stokesian  type  series  are 
nonuniformly  convergent,  and  are  only  valid  in  deep  and 
intermediate  water  depth. 

Boussinesg  (1877)  and  Korteweg  and  de  Vries  (1395) 
developed  the  cnoidal  wave  theories,  which  are  based  upon 
power  series  in  terms  of  wave  height  relative  to  water 
depth.  These  power  series  are  uniformly  convergent  in 
shallow  water,  but  are  invalid  in  deep  water. 

Baking  the  connection  between  the  two  above  theories, 
Goda  (1983)  developed  an  empirical  parameter  to  describe  the 
phenomena  of  wave  ncnlinearities  with  good  results.  The 
parameter  bridges  the  wave  steepness  in  deep  water  and  the 
Ursell's  parameter  Or  =  (a/d)  /  (d/L)  *  (with  a  the  wave 
amplitude,  d  the  water  depth,  and  L  the  local  wave  length) 


in  very  shallow  water,  spanning  the  fall  range  of  water 
waves. 

When  confronted  with  experimental  results,  Le  Hehaute  et 
al.  (1S68)  demonstrated  that  the  nonlinear  theories  were  not 
such  tetter  than  the  linear  wave  theory;  nonetheless  there 
is  a  need  to  clarify  the  extent  of  wave  nonlinearity  in  the 
sea  so  that  engineering  problems  can  be  solved  with  much 
more  confidence  and  accuracy. 

During  the  last  two  decades,  numerical  approaches  based 
on  truncated  Fourier  expansions  have  been  proposed  and  well 
verified  experimentally.  The  first  such  approach  was  devel¬ 
oped  by  Chappelear  (1961)  involving  the  use  of  the  velocity 
potential.  Dean  (1965)  used  the  stream  function  to  develop 
his  numerical  wave  theory,  which  was  computationally  simpler 
than  Chappelear* s  tcchnigue.  Dean's  stream  function  form 
satisfies  the  Laplace  equation,  the  Jcinematic  free  surface 
boundary  condition,  and  the  bottom  boundary  condition;  the 
parameters  in  the  stream  function  expression  are  chosen  by  a 
best  fit  to  the  dynaxic  free  surface  boundary  condition.  Von 
Schwinfi  and  Reid  (1972)  developed  characteristic  solutions 
to  finite  amplitude  waves  using  the  velocity  potential. 
Cokelet  (1977)  extended  the  method  originally  developed  by 
Schwartz  (1974)  to  allow  a  very  accurate  calculation  of  the 
characteristics  of  water  waves.  The  procedure  involves 
expressing  the  complex  potential  solution  in  a  Fourier 
series  and  represents  the  Fourier  coefficients  in  terms  of  a 
perturbation  parameter.  Rienecker  and  Fenton  (1980)  used  a 
finite  Fourier  series  to  describe  waves  by  solving  a  set  of 
nonlinear  equations  using  Newton's  method.  Recently,  le 
Hehaute  et  al.  (1984)  developed  a  parameterized  solution 
which  is  valid  for  all  practical  ranges  of  values  of  wave 
amplitude,  frequency  and  water  depth.  In  view  of  its  rela¬ 
tive  simplicity,  and  experimental  validity,  the  parameter¬ 
ized  solution  is  recommended  for  engineering  applications. 


The  objective  of  this  work  is  to  utilize  the  spectral 
model  by  le  Hehaute  et  al.  (1984)  to  coapare  with  wave  data 
acquired  at  Santa  Barbara  during  the  2-6  February,  1980. 
The  range  of  application  of  the  nodel  will  be  checked 
against  aeasured  wave  spectral  coaponents  and  individual 
waves.  The  wave  nonlinearity  and  Ur sell  nuaber  are  analyzed 
using  the  aodel.  Ccaaents  and  conclusions  are  nade  for 
further  study. 


1.  PABAHETEBIZED  SOIUTIOH  TO  THE  HOHLIHEAB  BATE  PBOBIEM 


1 .  Overview 

The  higher-order  Stokian  wave  theories  becoae  alge¬ 
braically  very  complex.  For  practical  purpose,  it  was 
desirable  to  develop  nave  theories  that  could  be  coaputed  to 
any  order.  Le  Hehaute  et  al.  (1984)  developed  a  parameter- 
ized  solution  to  the  aonochroaatic,  irrotational,  nonlinear 
vater  wave  over  a  horizontal  bottoa.  The  formulation  is 
expressed  in  a  closed  fora  in  terms  of  a  truncated  Fourier 
series.  All  coefficients  are  expressed  in  terns  of  simple 
analytical  functions  which  contain  three  paraneters,  namely, 
wave  height,  H,  wave  period,  T,  and  vater  depth,  d.  The 
number  of  terns  in  the  Fourier  series  is  proposed  paramet¬ 
rically  for  a  given  level  of  accuracy.  The  solution  is 
assumed  to  be  applicable  for  all  practical  ranges  from  deep 
to  shallow  water  waves,  i.e.,  d/Lo  >  0.005.  The  deep  water 
soluticn  for  d/lo  *  0.5  is  assumed  to  be  valid  for  d/Lo  > 
0.5. 


2.  Theoretical  Formulation 

As  in  the  periodic  water  wave  problem  the  formula¬ 
tion  has  linear  and  nonlinear  boundary  conditions  and  a 
linear  governing  differential  eguation. 

Choosing  a  coordinate  system  moving  with  wave  speed 
c  *  O'/  k  so  that  the  time  dependency  is  removed,  i.e.,  a 
staticnary  problem  ,  it  is  possible  to  assign  a  constant 
value  to  the  stream  function  at  the  free  surface.  The  vari¬ 
able  (J  is  the  wave  angular  f  rguency  (  O’  *  2  7T /T)  and  k  is 
the  wavenumber  (k  =  2  7T  /  L) .  The  solution  is  a  velocity 
potential  function  of  the  fora 


J*V< 


<£>  { (x-ct)  ,z]  or  <p[Q»z)  which  satisfies  the  Laplace 
equation  ; 


V20 


(2.1) 


The  variable  Q  =  kx  -  g"  t  denotes  phase,  and  the  z  axis  is 
vertical  positive  upwards  from  the  S.H.L.  (still  water 
level)  and  the  horizontal  x  axis  is  positive  in  the  direc¬ 
tion  cf  wave  propagation. 

The  kinematic  bottom  boundary  condition  is  written 


<f>z  •  •  0 

z*-d 


(2.2) 


The  dynamic  free  surface  boundary  condition,  is 
derived  using  Bernoulli* s  equation  written  in  the  dimen¬ 
sionless  form  at  z  =  T) 


i 


(2.3) 


where  H  is  the  wave  height,  TJ  is  the  water  wave  surface 
elevation  and  Be  =  as  arbitrary  Bernoulli  constant. 

The  kinematic  free  surface  boundary  condition  is 


vr- 


=  ’2(9) 


(2.4) 


The  e's  are  the  allowable  errors  which  should  be  small. 

The  measure  cf  how  well  these  two  boundary  condi¬ 
tions  are  satisfied  is  defined  by  E  t  ,  and  »  which  are 
the  mean  squared  errors  of  the  dynamic  and  kinematic  free 
surface  boundary  conditions  : 


V»v**k 


in  which  j  *  sample  at  evenly  spaced  phase  angles  (e.g. , 
with  5  °  interval)  along  the  wave  profile.  The  Bernoulli 
constant.  Be,  is  determined  such  that  the  variation  cf  the 
phase  speed  c  is  ainixized  over  a  wavelength.  For  an  exact 
solution,  2/  and  2%  would  be  zero. 

1  solution  in  the  fora  of  a  Halted  series  (or  trun¬ 
cated  Fourier  expansion)  is  assuaed 


<f>  a 

a  f 


coshlnk  (d  +z)]sin  nQ 


(2.7) 


in  which  <£)  *  diaensionless  potential  function;  An  *  a 
derived  coefficient;  n  is  an  integer  and  N  =  the  number  of 
terms.  All  the  internal  wave  field  terns  (velocity  compo¬ 
nents,  dynamic  pressure  and  acceleration  components)  can  be 
directly  determined  from  this  equation.  Eg.  2.7  satisfies 
both  the  Laplace  eguation  (2.  1)  and  the  kinematic  bottom 
boundary  condition  (2.2)  exactly. 

The  eguation  for  the  free  surface  correponding  to 
eg.  2.7  is: 


N 

V  *  S  *  An  sinhlnk  rd +  ri)]cos  n  f) 

2a  20  2<J*/_ _  L  J  (2.8) 

ns  i 


in  which  K  =  an  integration  constant  satisfying 


At  the  wave  crest  {  Q  *  0)  ,  7f  *  7)e  and  the  vave  trough 

(  8  *  J  ,  77  *  Tjt  -Therefore 


Ik  +  Lad  =  /  + 

2  e 


(2.  10) 


in  which  2  ^  is  to  be  determined  and  aust  be  snail. 

le  Hehaute  et  al  (1984)  defines  the  phase  velocity  c 
or  vavenunber  X.  in  accordance  with  the  Hamiltonian  varia¬ 
tional  principle  of  xiniaun  energy,  which,  in  practice,  is 
obtained  through  minimizing  2 1  with  respect  to  the  waven- 
unber,  k 


<9  £, 

d  *  "  0  <2-"> 


Such  a  system  of  equations  can  eventually  be  solved  as  a 
function  of  phase  velocity. 

The  solution  to  the  problem  consists  in  providing  an 
analytical  solution  which  describes  all  the  wave  character¬ 
istics  as  functions  of  only  three  parameters  :  wave  height, 
H,  wave  period,  T,  and  water  depth, d.  These  paraneters  are 
grouped  into  two  dimensionless  parameters 


Depth  parameter: 


2  ITS 


(2.12) 


Have  steepness  parameter: 


(2.  13) 


is  which  lo  is  the  linear  deep  water  wavelength  and  a  *  H/2 
is  the  wave  aaplitude. 

Having  specified  the  initial  boundary  value  prcblea 
and  analytical  foras  cf  egs.  2.3  and  2.4,  the  problea  now 
consists  of  expressing  all  unknowns  ;  An,  k,  K,  and  N  (eg. 
2.8),  in  teras  of  the  two  paraaeters  <5  aad  y  ,  such  that 
a  desired  level  of  accuracy  is  achieved.  These  unknowns  are 
deter ained  by  satisfying  the  boundary  criteria  as  closely  as 
possible.  2g.  2.11,  which  specifies  the  phase  velocity  and 

wave  nuaber,  must  also  be  satisfied.  The  aethod  of  solution 
is  one  cf  trial  and  error  involving  aultiple  iterative 
processes  as  expressed  in  the  flow  chart  in  Figure  1. 

For  practical  purposes,  the  final  foraulation  nay 
be  presented  in  the  following: 

INPUT  EATA  :  Given  E,  T,  and  d  then 

7Ty  **=27r/t*  ,  *=27r/t  ^  o-  =  27r/r 

and  ~Y  ss H /  L  ai ,  ~.°-995 

1  T  mnlLo  ,Als0j  p^y /tonh  (4.5888  ) 

and  (T2  =  j  k*  tanh  k*d  . 

(The  syntol  *  denotes  the  linear  solution) . 

FORM  CF  SGLOTIONS: 

N 

$=  (4 +  l^sin„e  (2.14) 


T1  =  K-  ^-"*['>*(^  +  77)]  CO*  (2  15) 
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INTEGRATION  CONSTANT:  For  application  purposes,  a  parame¬ 
terized  form  of  K  is  determined  by 


K  =  -  1.764  j  l £  p“  SS§e'77 


<2.  16) 


The  use  of  K  as  expressed  above  generally  yields  larger 
values  cf  E*s  than  the  value  of  K  obtained  by  eg.  2.9 
directly. 

NONLINEAR  WAVE  NONBEE  S  k  *  2  TT/L 
1.  0.005  <5  <0.1 


*=»*'= _ i 

r  .  nl.097  _  _ 


[  /  +  l.879p,^^0.094p^Jn(2,lo28)] 


2.  0.1  <  5  <  0.5 


k  as  t*a 


/  ,  3.466  nl%$  rs0.3 15  1 

tanh  k*d  $  ®  J 


(2.  17) 


(2.  18) 


The  solution  valid  for  Q  =  0.5  is  also  considered  valid  for 
5a“"  0.5  provided  that  all  the  parameterized  functions  which 
follow  are  determined  with  5=0-5  irrespective  of  their 
real  value;  this  means  that  deep  water  solutions  corre¬ 
sponding  to  5  >  0.5  are  not  influenced  by  water  depth. 

COEFFICIENTS: 

Coefficient^  A,: 


A,  =  A*(/  -  R  ) 


(2.  19) 


cosh  k*d 

1.  0.005  «  S  <0.1 


s  .  «;  (S.sszp )*I 

in  which  a«y  *  0.8132  exp  {-0.0153  Jin  (1  O'3 
b}  *  0.0747  exp  {1.  1057  Jin  <10J  5 
2.  0.1  <  5  ^  0.5 

R  =  0.3  +  (-0.3  +  a” )  explb*  [ln  ( 10  5)  J  /  , 


in  which  a*J  *  17.30^ 


2.38 


b*J  *  1.389  tanh,(24.0S^}/*8'36  ) 

n3.04 

c'j  *  1.3676  exp  <-179.0/3  ) 

Coefficient  A,  : 

1.  0.005  ^  5  <  0- 1 


A2  =  0.497  A^xp^fS.BSZf]  )b2  j 

in  which  a^  =  -0.00933  exp{1.494  J  In  ( 1 03  i 
b*2  -  -0.953  exp  {-0.029  [lnflO^S 
2.  0.1  <  5  4  0.5 

M 

A2  *  2.5x10  5 Af  txp[a2exp(b2  S^)] 


in  which  a "  »  10.547  exp  (-0.0 134 ft 


-o.9g 


) 


b»  «  —  1 6 ♦  1 5»  318  exp(-Q.0017^J 
c«  *  5-4.328  exp  (-0.01  M/]'104  ) 
Coefficient  Aj  : 

1.  0.305  «  S  4  0-  1 


A3  rn  0.672  A2  exp  £  a'3  (5.8*2/])**] 

7& 

in  which  a*j  *  -0.021  exp  {1.355  |ln(103  6  )J  °  ) 

b'3«  -0.859  exp {-0.0747  [ln(1036  )^'Z*  } 
2.  0.1  <64.  0.5 


A3  *  A2  **P  (o') 


(2.26) 


ro.52  **-o.  41 8 

in  which  a"  *  -3.5  37  0  (J 

w 

Coefficient  An  (n  >  3) 

ro.S2  _oA,8 

An  -  An_,  exp  (-3.5370  f]  ) 


(2.27) 


3.  ISalxtjgal  lalidjti 

The  validity  cf  this  theory  is  defined  through  the  E 
values  as  derived  by  egs.  2.3,  2.4,  2.5,  2.6,  2. 1C  and  2.11 
by  order  of  importance.  The  magnitude  of  the  E*s  depends  on 
the  nucber  of  terns  H  used  in  the  series.  3eyond  a  certain 
value  of  N,  it  is  found  that  the  E  values  tend  towards 
asymptotic  values,  indicating  that  they  are  not  reduced  by 
adding  more  terns,  i.e. ,  the  series  is  bounded.  With  this 


property  in  the  solutions,  the  nuaber  of  teras  >  reguired 
can  be  expressed  in  parameterized  fora  by  : 


II  * 

integer  part  of(a n2(5.38 

P  *  °-2)} 

(2.  28) 

n/* 

1  ♦  19.13  exp (-0.344  |ln 1 10 3 

S  )]C1.‘2>} 

(2.29) 

nz~ 

32.72  exp (-0.212 |ln (10*  5 

)]<».*»>} 

(2.30) 

le  Mehaute  et  al,  (1964) 

showed  that  the 

parameter 

ized  solnticn  approxiaates  the  free  surface  dynamic  boundary 
condition  as  well  as,  if  not  better,  than  any  other  existing 
theory.  It  also  approximates  the  kinematic  condition  gener¬ 
ally  better,  except  for  the  stream  function  theory  devel¬ 
oped  by  Dean  (1974)  by  virtue  of  its  definition. 

Truncated  Fourier  series  do  not  represent  limit 
waves  well  since  the  cusp  at  the  wave  crest  reguires  an 
infinite  number  of  terms  as  was  demonstrated  by  Dean  (1974). 
Therefore,  the  accuracy  of  the  parametric  solution  decreases 
rapidly  as  the  wave  heights  approach  limit  wave  conditions, 
which  is  where  Ur  is  large. 

B.  MIASIJBE  OF  NOHLIBf ABITT 

The  parameter  which  represents  the  wave  nonlinearity  in 
deep  water  is  the  wave  steepness,  or  a  measure  of  the 
surface  elevation  slope.  The  most  common  representation  is 
either  the  ratio  of  wave  height  to  wavelength  H/L,  or  the 
parameter  aX. 

In  the  region  of  shallow  water  waves,  the  Stokian  wave 
profile  for  higher  orders  is  an  expansion  using  the  nondi- 
mensicnal  term 


Ur  *  (a/d)/(d2/L2) 


(2.31) 


as  the  perturbation  parameter.  Ibis  tera  is  called  the 
Orsell  parameter  and  has  been  shown  to  govern  the  transform 
nation  of  water  waves  in  shallow  water.  For  exaaple,  it  has 
been  deaonst rated  that  the  nonlinear  shoaling  of  water  waves 
is  predicted  as  the  function  of  the  Orsell's  paraaeter. 

When  Or  «  1,  the  linear  snail  aaplitude  wave  theory 
applies.  In  principle,  nore  and  aore  terns  of  the  power 
series  are  required  in  order  to  keep  the  sane  relative  accu¬ 
racy  as  the  Orsell  paraaeter  increases. 

Ihe  Orsell  paraaeter  is  a  usefull  guide,  but  is  not 
necessarily  sufficient  for  judging  the  relative  importance 
of  the  nonlinear  effects.  &  qualitative  idea  of  the  iwpcr- 
tance  of  nonlinearity  is  also  given  by  the  ratio  of  spectral 
values  of  each  harmcnic  frequency  to  the  spectral  value  of 
the  fundamental  frequency.  However,  the  nonlinear  effects  of 
shoaling  transformation  are  cumulative,  so  that  the  ncnli- 
nearities  may  be  locally  weak,  but  significant  cross- 
spectral  transfer  can  occur  if  the  shoaling  region  is  wide 
in  comparison  to  the  nonlinear  interaction  distance. 


III.  EXPBBIHEHTAX  PBOCEDOBBS 


Dave  and  velocity  data  acquired  at  Santa  Barbara, 
California,  during  tbe  2nd,  3rd,  4th,  5th  and  6th  February 
1980,  as  part  of  the  Nearshore  Sediment  Transport  Study,  are 
used  to  be  coapare  with  the  paraaeterized  aodel.  The  field 
data,  considered  in  this  work,  include  only  measurements 
aade  by  the  current  meters  outside  the  breaking  zone. 
Current  meters  C07X,  C03X  and  C91I  (see  Figure  2)  are  used 
on  the  2nd  and  6th  of  February;  current  meters  C03X,  C01X 
and  CODX  are  used  on  the  3rd  of  February,  and  CO IX  and  CODX 
are  used  on  the  4th  of  February. 

For  each  day  and  each  instrument  the  following  statis¬ 
tics  were  calculated:  average  wave  height,  F,  the  maximum 


wave  height,  Haax,  tbe  significant  wave  height,  Hl/3,  the 
root  mean  squared  wave  height,  Hrms,  and  the  average  of  the 
heights  cf  the  1/10  highest  waves,  Hi/10,  as  well  as  the 
associated  water  depths,  and  the  surface 
elevation  spectrum. 


A.  HAVE  SPECTBAL  A1A1YSIS 

The  surface  elevation  spectra  were  calculated  from  63 
minute  current  meter  records.  The  water  particle  velocity 
components  measured  by  the  current  meters  were  recorded  on  a 
special  receiver/tape  recorder,  digitally  low-pass  filtered 
and  sampled  at  2  samples/sec.  The  data  were  then  high-pass 
filtered  at  0.05  hz  to  exclude  low  frequency  variations. 
These  current  data  were  used  to  infer  wave  heights.  The 
velocity  signals  were  convolved  using  linear  wave  theory  to 
obtain  surface  elevations.  The  complex  Fourier  spectra  of 
the  horizontal  velocity  components  were  first  calculated  and 


vecrt  orially  added  to  give  7  (£)  .  The  complex  surface  eleva¬ 
tion  spectrun,  2(f),  vas  calculated  applying  the  linear  wave 
theory  transfer  function  H(f) 

X(f)  =  H  {£)  *  7(f)  (3.1) 
where 

H  (f )  =  (sinh  kd)/{(Jcosh  k(d+z)}  (3.2) 

with  z  =  measurement  cf  elevation. 

The  calculated  energy  density  spectra  are  to  be  compared 
with  the  parametric  model.  Since  the  parametric  model 
defines  all  the  truncated  Fourier  series  coefficients  as 
simple  analytical  functions  in  terms  of  the  wave  height  H, 
wave  period  T,  and  water  depth  d,  it  was  necessary  to  trans¬ 
form  the  energy  density  spectra  of  the  sea-surface  elevation 
into  a  discrete  (spikes)  spectra  as  shown  in  Figure  3. 

The  discrete  spectra  were  calculated  by  first  finding 
the  fundamental  freguency  defined  as  the  frequency  associ¬ 
ated  with  the  highest  energy  density  peak.  The  harmonics 
were  defined  as  the  integer  mutiples  of  this  fundamental 
frequency.  The  bandwidths  associated  with  the  energy  peaks 
were  defined  at  the  frequencies  corresponding  to  the  half 
power  points  of  the  energy  density  values.  The  variance  of 
the  peaks,  calculated  as  the  area  of  the  energy  density 
spectrum  between  half  power  points,  describe  the  height  of 
the  spectral  spike.  The  analysis  was  carried  out  to  the 
highest  harmonic  for  which  half  power  points  could  be 
defined.  The  reason  for  defining  the  bandwidths  at  the  half 
power  points  is  because  the  half  power  values  correspond  to 
the  level  at  which  the  spectral  components  are  independent 
of  one  another.  It  was  noted  that  the  bandwidth  increases 
approximately  with  the  number  of  harmonic. 


Assuming  Tj  3  0 ,  the  varxance  of  the  wave  surface  xs 

given  by 


var  * 


f/  »2-' 

JT 


(3.3) 


The  parameterized  variance,  or  energy,  can  be  defined  fron 
equation  2.15  in  the  discrete  aode 


(3.4) 


with  n  =  1,2,3...B. 

Ihe  wave  height  (and  wave  aaplitude)  used  to  calculate 
the  parametric  energy  line  spectrum  was  the  average  wave 
height,  IT, because  it  was  the  wave  height  parameter  which 
offered  best  fitting  to  the  field  data. 

Ccsparisons  for  all  days  were  made  between  the  converted 
actual  line  spectrua  and  the  parametric  line  spectrum  calcu¬ 
lated  using  the  input  parameters  of  average  wave  height  H, 
average  water  depth  d,  and  the  period  corresponding  to  the 
defined  fundamental  frequency  (1  =  1/f) ,  calculated  for  the 
68  minute  measurements. 


B.  ISEIVIDUAL  RAVE  A1ALYSIS 

A  second  approach  was  to  compare  individual  waves  meas¬ 
ured  in  the  field  with  the  parametric  model  predictions. 
Ihe  filtered  and  digitized  velocity  signals  were  convolved 
using  linear  theory  tc  obtain  surface  elevations.  The  indi¬ 
vidual  waves  were  defined  using  the  zero-up  crossing  method 
as  explained  in  Figure  10.  The  highest  maximum  and  the 
lowest  minimum  of  the  surface  elevation  within  a  period 
interval  define  the  crest  and  the  trough  of  a  wave.  A  wave 
height  Hi  is  defined  as  the  range  of  Tj  (t)  in  that 
interval,  the  period  is  the  time  interval  between  two 
consecutive  zero-up  crossings  cf  Tj  (t)  . 


A  Fourier  analysis  was  perforaed  on  each  individual  wave 
using  a  Fast  Fourier  Trans fora  algorithm.  The  parametric 
coefficients  associated  with  the  surface  elevation  formula 
were  also  determined  for  each  wave  using  the  measured  Hi, 
Ti,  and  di.  Statistics  of  Fourier  coefficients  from  field 
data  and  the  parametric  model  were  generated  for  each  day 


w 

Tn  »  (1/Hw)  *  ^An(i) 


(3.5) 


(where  Nw  =  number  of  analyzed  waves).  For  each  run  compar¬ 
isons  were  made  between  this  approach  and  the  wave  spectral 
analysis.  Since  the  average  wave  period  was  about  15 
seconds,  the  total  nueber  of  waves  in  each  68  minutes  record 
was  about  300,  which  give  reasonable  wave  statistics. 


C.  S0HHAB1 

Both  analyses  were  based  on  the  similarity  between  the 
Fourier  series  representation  of  the  sea  surface  elevation 


TJ  (t)  *  aO  ♦  an  cos  nQ 


(3.6) 


(with  aO  »  mean)  and  the  sea  surface  elevation  described  by 
the  parametric  model  as  defined  by  the  equation  2.15.  The 
coefficients  to  be  compared  are  the  an*s  for  the  Fourier 
series  defined  by  the  Fourier  integrals  and  the  Le  Nehaute 
coefficients 

-  (1/  CT  **2)  ♦a*g*h*  An  *  sinh  (nk  (d+77  )} 
where  An  are  defined  parametrically  as  function  of  the  depth 
parameter  (5  and  wave  steepness  y  ,and  TJ  defined  itera¬ 
tively  from  equation  2.15. 


The  spectral  data  analysis  is  a  linear  approach.  The 
spectral  coaponents  are  assaaed  independent  of  ?ach  other. 
The  analysis  considers  that  the  spectral  coaponent  phase 
angles  at  each  frequency  are  randoa  (putting  the  preblea  in 
a  probabilistic  or  stcchastic  setting).  On  the  other  hand, 
the  individual  wave  analysis  is  also  a  linear  approach,  but 
the  phase  relationship  between  haraonics  is  presuaed  deter-* 
ainistic  (i.e.  phase  coupled  with  the  fundaaental  frequency, 
as  evidenced  by  the  observed  peaked  crests  and  flattened 
troughs) . 


IT.  HE SUITS 


1.  SPECTRAL  AI1LTSIS 

1 .  lonlinear  Statistics 

A  typical  energy  density  spectrum  for  all  analyzed 
instruments  is  presented  in  Figure  2  (upper  panel) .  The 
presence  of  strong  haraonics  in  the  spectrua  indicates  the 
iaportance  of  nonlinearities  of  the  waves  in  this  region  of 
shallow  water.  A  qualitative  idea  of  the  iaportance  of 
nonlinearity  is  given  by  the  ratio  of  the  harnonic  energies 
to  the  priaary  frequency  energy  E(nf)/E(f),  where  n  =  2,  3 
...  are  the  harnonic  frequencies.  Due  to  the  difficulty  of 
defining  the  half  power  points  for  the  haraonics  higher  than 
the  third,  only  the  ratios  involving  the  first  and  second 
harmonics  are  presented. 

Tables  I  -  V  give  the  ratios  for  all  five  days 
studied  for  both  the  paraaeterized  aodel  (denoted  L. n.)  and 
the  field  data  (denoted  by  F.).The  water  depth  and  the 
Ursell  parameter  associated  with  each  measurement  are  also 
presented.  The  range  of  water  depths  is  from  about  670  cm 
to  138  ca.  The  range  cf  the  Ursell  parameter  is  from  abcut  6 
to  12C. 

For  large  Ursell  numbers  (  Ur  >  75  ) ,  the  parameter¬ 
ized  acdel  overestimates  the  iiportance  of  the  energy  asso¬ 
ciated  with  the  harmonics  relatively  to  the  fundamental 
frequency  energy.  The  same  ratios  for  the  field  data  are  20 
-  60  percent  lover.  For  small  Ursell  numbers  (  Or  <  25  ) , 
the  parameterized  aodel  underestimates  the  importance  of  the 
harmonics  relatively  to  the  fundamental  frequency  energy. 
The  same  ratios  for  the  field  data  are  0.5  -  40  percent 


The  best  fitting  between  the  parameterized  aodel  and 
the  field  data  occurs  at  intermediate  Orsell  numbers  (  25  < 
Or  <  75  ) .  It  is  found  that  the  parameterized  aodel  energy 
ratios  compared  with  the  field  data  energy  ratios  have  an 
error  of  less  than  20  percent  for  the  cases  of  intermediate 
Orsell  number.  Some  exceptions  were  observed  to  the  above 
indicated  percentages,  but  the  energy  contribution  at  these 
frequencies  is  very  small. 

2.  local  Spectral  Comparisons 

In  absolute  comparison  was  made  between  the  parame¬ 
terized  aodel  spectral  energy  peaks  and  field  data  spectral 
energy  peaks  as  a  function  of  frequency.  Typical  results  for 
three  different  Orsell  numbers  are  presented  in  Figures  4 
and  5.  In  regions  of  strong  nonlinear  effects  (  large 
Orsell  number)  ,  the  parameterized  model  underestimates  the 
energy  at  the  fundamental  frequency.  The  energy  drops 
linearly  from  the  fundamental  frequency  to  the  harmonics,  in 
the  parameterized  model,  while  field  data  spectra  fall  cf£ 
much  faster.  In  regions  of  small  Orsell  numbers  the  parame¬ 
terized  model  concentrates  almost  all  the  energy  in  the 
fundamental  frquency,  overestimating  it  relative  to  the 
field  data.  The  best  agreement  between  the  parameterized 
spectral  energy  peaks  and  field  data  energy  peaks  are 
observed  for  intermediate  Orsell  numbers. 

The  Fourier  coefficients  used  to  evaluate  field  data 
spectra  are  plotted  against  the  Le  Nehaute  spectral  coeffi¬ 
cients  (  1  /  (J  **2)  agkA  sinh[nk(d  ♦77)}  (Figs.  6-8).  A 

reasonable  correlation  is  found  for  the  first  three 
coefficients. 

The  spectral  analysis  approach  is  summarized  by 
plotting  the  ratio  of  the  Le  Mehaute  spectral  peaks  to  field 
data  spectral  peaks  against  the  (Jrsell  parameter  (Fig.  9) . 
The  ratio  for  the  fundamental  frequency  energy  decreases  as 
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the  Orsell  parameter  increases  emphasizing  hov  the  parame¬ 
terized  model  underestimates  the  fundamental  frequency  spec¬ 
tral  peak  in  situations  of  strong  nonlinearities  (  large 
Orsell  parameter).  Ihe  ratios  for  both  the  first  and  second 
harmonic  frequency  energies  increase  as  the  Orsell  parameter 
increases;  it  must  be  noted,  as  indicated  earlier,  that  the 
data  cnly  ranges  in  Orsell  numbers  from  about  6  to  120. 

B.  IISrVIDOlL  HITS  11A1ISIS 

Individual  waves  mere  identified  using  the  zero-up-cross 
method.  A  Fourier  analysis  was  performed  on  each  individual 
wave.  Ihe  parameterized  spectral  coefficients  were  also 
determined  for  each  individual  wave. 

Ihe  Fourier  coefficients  associated  with  the  fundamental 
frequency,  first  and  second  harmonics  are  plotted  separately 
against  the  corresponding  parameterized  coefficients.  Ihe 
coefficients  at  the  fundamental  frequency  (Fig.  14)  are 
reasonably  well  correlated  (correlation  coefficient  - 
0.80).  Ihe  Fourier  coefficients  plotted  against  the  Le 
Hehaute  coefficients  follow  a  straight  line  with  intercept 
at  the  origin  and  slcpe  of  about  one.  The  group  of  plotted 
points  with  a  greater  slope  (approximately  2)  are  associated 
with  measurements  in  shallower  waters  with  higher  Orsell 
numbers  (see  Figs.  12-16).  The  graphs  of  the  the  first  ar.d 
second  harmonic  coefficients  (Figs.  17-18)  show  similar 
characteristics  as  the  fundamental  frequency;  however,  the 
correlation  between  the  coefficients  is  poorer. 

Ihe  ratio  of  the  le  Hehaute  coefficients  to  the  Fourier 
coefficients  is  also  plotted  against  the  Orsell  parameter 
Ihe  ratio  for  the  fundamental  frequency  coefficients 
converge  to  a  bound  value  of  0.5  at  large  Orsell  numbers 
(Fig.  19).  For  Orsell  numbers  less  than  75  the  ratio  of  the 


le  Hehaute  to  Fourier  coefficients  is  about  1  with  the 


exception  of  a  few  points  near  Ur  *  0  which  were  identified 
to  he  related  with  noisy  signals.  For  large  Ursell  numbers 
the  Fourier  coefficients  are  about  twice  the  Le  Hehaute 
coefficients,  again  showing  the  tendency  that  the  parameter- 
ized  aodel  has  to  underestiaate  the  fundaaental  peak  energy 
in  regions  of  strong  nonlinearities  (large  Ursell  nuabers). 

The  ratios  for  the  first  and  second  haraonic  coeffi¬ 
cients  (Figs.  19  and  20)  show  a  bounded  value  near  unity  as 
the  ursell  paraaeter  increases.  Hoisy  signals  near  Ur  =  0 
are  also  observed. 

C.  CCBPABISOHS  BETIZZI  1PPR0ACHBS 

Both  linear  (spectral  analysis)  and  nonlinear  (indi¬ 
vidual  wave  analysis)  approaches  lead  to  siailar  results. 
The  paraneterized  acdel  is  observed  to  overestiaate  the 
fundaaental  frequency  energy  in  regions  of  snail  Ursell 
nunber  (weak  nonlinearities)  and  underestiaate  it  in  regions 
of  large  Ursell  nunber.  The  individual  wave  analysis  pres¬ 
ents  a  liait  tendency  for  large  Ursell  nuabers,  with  the 
fundaaental  frequency  Fourier  coefficients  about  twice  those 
of  the  paraneterized  spectral  coefficients.  The  better 
correlation  observed  between  the  Fourier  coefficients  and 
the  le  Hehaute  coefficients  in  the  individual  wave  analysis 
as  conpazed  with  spectral  analysis  approach  is  due,  perhaps, 
to  the  definition  of  half-power  criteria  or  to  the  averaging 
processes  involved  in  calculating  the  wave  spectra. 

The  reason  for  the  sinilarity  of  results  for  the  two 
approaches  can  be  seen  by  coaparing  the  spectra  calculated 
in  the  usual  manner  with  a  "nonlinear  spectrum".  &  nonli¬ 
near  spectrum  was  calculated  using  the  Fourier  amplitude 
coefficients  calculated  for  the  individual  waves.  Using  the 
same  handwidths  as  the  ordinary  spectrum,  amplitude  coeffi¬ 
cients  were  sorted  into  the  respective  frequency  bands  and  a 
nonlineaz  spectrum  calculated  as 


Gal  (f)  *  (1/*w)  ^  1?  (f)  /  2  4f 


(4.1) 


■here  the  coefficients  Ai(f)  are  centered  on  frequency  f 
within  bandwidth  A  f.  The  total  nuaber  of  waves  is  denoted 
by  Nw.  in  ezaaple  of  the  coaparison  of  the  nonlinear  spec- 
trua  calculated  in  this  Banner  with  the  linear  spectral 
aethod  is  shown  in  Figure  21.  The  two  spectra  are  siailar 
and  bcth  show  strong  haraonic  peaks.  The  spectral  energy  in 
the  ordinary  spectrua  tends  to  be  saeared  aore,  whereas  the 
peaks  and  valleys  in  the  nonlinear  spectrua  are  aore  sharply 
defined. 


The  coefficients  of  the  parameterized  solution  tc  the 
monochromatic,  irrotational,  nonlinear  wave  over  a  hori¬ 
zontal  bottom,  as  developed  by  Le  Behaute  et  al.  (1984)  ,  are 
expressed  in  teras  cf  simple  analytical  functions  of  wave 
height,  H,  uave  period,  T,and  vater  depth,  d.  All  quanti¬ 
ties  explicitly  derived  froa  the  potential  function 
(velocity  coaponents,  dynaaic  pressure,  acceleration,  i.e. , 
the  internal  Have  field)  can  be  directly  deterained.  The 
relative  simplicity  of  the  calculations  of  the  parameterized 
solution  recommends  it  for  engineering  applications. 

The  Have  spectra  determined  from  the  parameterized  solu¬ 
tion  shou  good  agreeaent  with  the  measured  Have  spectra  from 
field  data  for  conditions  of  intermediate  Orsell  numbers  (25 
<  Or  <  75) .  For  saall  Orsell  numbers,  the  parameterized 
solution  concentrates  almost  all  the  energy  in  the  funda¬ 
mental  freguency  representing  a  sinusoidal  vaveform.  Under 
these  conditions,  the  field  vave  spectra  always  contained 
relatively  less  energy  associated  with  the  funda cental 
freguency  and  more  energy  distributed  to  the  harmonics. 

For  large  Orsell  numbers,  the  parameterized  solution 
describes  a  waveform  with  more  peaked  crests  and  flatter 
throughs  relative  to  the  measured  waves,  since  it  underesti¬ 
mates  the  energy  associated  with  the  fundamental  freguency 
and  overestimates  the  energy  associated  with  the  higher 
harmonics. 

The  ratios  between  the  Le  Hehaute  and  field  data  spec¬ 
tral  coefficients  associated  with  the  1st  and  2nd  harmonics 
present  a  bound  licit  tendency  of  1  for  large  Orsell 


numbers,  indicating  good  agreement  between  the  two  processes 
for  these  conditions.  The  sane  ratios  for  the  coefficients 
associated  with  the  the  fandaaental  frequency  show  a  hound 
limit  tendency  of  0.5  ,  indicating  that  the  first  coeffi¬ 
cient  of  the  paraaeterized  model  needs  to  be  corrected  for 
strong  (Or  >  75)  nonlinear  effects. 

Further  coapariscns  should  be  aade  for  different  sea 
conditions  to  deteraine  the  corrections  to  the  parameterized 
model  coefficients.  However,  it  must  be  noted  that  a  lini- 
tation  of  the  parameterized  solution  is  that  its  accuracy 
decreases  rapidly  as  the  wave  hei^uts  approach  limiting  wave 
conditions  due  to  the  use  of  a  truncated  Fourier  series. 
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TABLE  II 

Ratio  of  Harmonic  Primary  Freguency  Energy 
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Figure  6  Fundamental  Frequency  Coefficients  Correlation 
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Spectral  Peaks  Comparisons  against  Orsell  N naber 


Figure  11  Coefficient  Ratio  at  Fundamental 
Frequency  against  the  Orsell  Munber. 
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Coefficient  Ratio  at  1st  Haraonic 
gainst  the  Orsell  Nuaber. 


Figure  13  Coefficient  Ratio  at  2nd  Haraonic 
against  the  Ursell  Nuaber. 


Individual  Waves 


1c2rrela^ion„o£  the  Coefficients  at 
Fundaaental  Frequency  for  all  the  Instruments 


Correlation  of  the  Fundamental  Frequency  Coefficients 
for  the  Current  Meter  CQ7  -2Fet.80. 


correlation  of  the  Fundamental  Frequency  Coefficients 


Figure  17  Correlation  of  the  Fundamental  Frequency  Coefficients 
for  the  Current  Meter  C01  -  2Feb. 80. 


Correlation  of  the  Fundamental  Frequency  Coefficients 


Figure  20  Correlation  of  the  Coefficients 
Haraonic  for  all  the  Instruaents. 
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